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Abstract We reduced the problem of constructing bi-Hamiltonian structure 
in three dimensions to the solution of a Riccati equation in moving coordinates 
of Serret-Frenet frame. We then show that either the linearly independent 
solutions of the corresponding second order equation or the normal vectors of 
the moving frame imply two compatible Poisson structures. 



1 Introduction 

The discovery of completely integrable nonlinear evolution equations as well as 
the algebraic and geometric structures associated with them has triggered an 
intensive search of finite dimensional dynamical systems resembling the similar 
properties. The bi-Hamiltonian structure as an underlying geometrical frame- 
work for complete integrability provoked the revival of Poisson structures of 
finite dimensional dynamical systems (see [1] and the references therein for de- 
tails and comparison). 

Several works [2]- [9] on construction of conserved quantities, on Hamiltonian 
structures and on integrability of three dimensional systems have led to a sys- 
tematic investigation using Poisson geometry, Frobenius integrability theorem 
and unfoldings of foliations [10], [11]. We presented correspondence between 
Poisson structures and integrable one forms and utilized this to obtain criteria 
for local and global existence of Poisson structures. We also obtained the local 
result that any two Poisson structures can be made into a compatible pair to 
form a bi-Hamiltonian structure. 

The restrictions on Poisson matrices imposed by the Jacobi identity and 
the compatibility condition for two of them to form a bi-Hamiltonian structure 
are the most serious conditions requiring deliberate actions against their simple 
presence as one scalar equation in three dimensions [12]- [15]. In [11], using 
an invariance property, we reduced the Jacobi identity to a nonlinear equation 
in ratios of components of the Poisson matrix. This scalar equation in one 
unknown function was shown to contain sufficient information for constructing 
the Poisson structure completely and was recognized to be the Riccati equation 
in [16]. 

The possibility to determine the Poisson matrix by a single function straight- 
ens out the difficulties in general Hamiltonian systems which does not fall into 
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the classes of canonical Hamiltonian or Lie-Poisson (i.e. with linear Poisson 
structure) equations, arising from the absence of coordinates similar to the 
canonical Darboux coordinates of symplectic geometry. The Darboux-Weinstein 
theorem [17] describes the local structinc of a Poisson manifold as a space fo- 
liated by symplectic submanifolds. The foliation depends on the rank of the 
structure which is an invariant of the Poisson matrix. As a result, the Poisson 
matrix consists of a constant submatrix whose rank is the rank of the Pois- 
son structure and some additional nonlinear part. In three dimensions, one 
must solve for at least one unknown function to determine the Poisson matrix 
completely. 

The role of the constants of motion, in particular the Casimirs, in lineariza- 
tion and integration of the equations associated with the Jacobi identity were 
also discussed in [12]-[16]. Yet, there is no direct relations of the Jacobi identity 
and the compatibility condition to the theory of linear differential equations 
which may be one of the elegant ways to avoid some deceptive conclusions from 
these simple looking differential equations ignoring their nonlinear character. 
The main source of these confusions is to endeavor to exploit local criteria for 
global results without questioning any suspicious obstructions such as the one 
we have found for the Darboux-Halphcn system [11]. 

In [11] we also observed through several examples that some coordinate 
transformations may cast the dynamical systems into a form where, in spite of 
much higher degree of nonlincarity, the integration for the conserved quantities 
and hence the manifestation of the bi-Hamiltonian structures become more ef- 
ficient. The non-covariance of the general Hamiltonian formulation is a source 
of the common belief that the existence of Hamiltonian structure and the intc- 
grability of dynamical systems rely heavily upon the coordinates in which they 
are represented as well as the prefered parametrization of the solution curves. 
Such coordinates were found to be important in numerical integrations as well. 
It is shown in [18] that, a necessary condition for numerical solution algorithms 
to preserve the conserved quantities of dynamical systems is their resemblance, 
as products of a skew symmetric matrix and a gradient vector, to Hamiltonian 
systems. See also [19] for numerical schemes applied to some concrete examples 
of systems under consideration. In three dimensions, the Nambu mechanics [2] 
is the only and generic (up to a conformal factor) framework which enables us 
to identify such coordinates. Namely, the Nambu structure is a manifestation 
of the bi-Hamiltonian structure, hence integrability, in a frame with coordinate 
vectors consisting of the dynamical vector field and gradients of two conserved 
Hamiltonians [11]. 

The Nambu representation obviously requires the integration of the system 
for Hamiltonian functions. In our study of the Darboux-Halphcn system, we 
were able to obtain obstraction for the global integration of such quantities [11]. 
To our knowledge, this is the only example which exhibits, along with a rich 
geometric structure, differentiation between local criteria and global availability. 
On the other hand, the local version of the Nambu mechanics has not yet been 
appeared in the literature. 

In this work, we shall show that the coordinates associated with the Serret- 
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Frcnct frame is the one we sought for three dimensional systems. In these 
coordinates the Jacobi identity hnearizes through the Riccati equation and the 
compatibiUty follows from the Hamilton's equations. Obstructions to the global 
constructions of the bi-Hamiltonian structures are encoded in the helicities and 
the cross-helicity of the unit vectors spanning the normal plane to the vector 
field associated with the dynamical system. We shall construct a bi-Hamiltonian 
moving frame for the local Nambu representation. 

In the next section we review the properties of bi-Hamiltonian systems in 
three dimensions. In section three we introduce the Serret-Prenet frame as- 
sociated with a vector in three dimension. We express the Jacobi identity in 
Serret-Frenet frame and show that in moving coordinates it reduces to a Riccati 
equation. In section four, we shall show that Poisson structures constructed 
from the solutions of the Riccati equations and/or the normal vectors arc all 
compatible via Hamilton's equations of motion. We then conclude the existence 
of bi-Hamiltonian structures. 



2 Hamiltonian Systems in Three Dimensions 

Following [1 1] , we shall summarize the necessary ingradients of the bi-Hamiltonian 

formalism in three dimensions. For x = {-t*} = (x, y, z) G M^, t G R and dot de- 
noting the derivative with respect to we consider the autonomous differential 
equations 

x = v(x) (1) 

associated with a three-dimensional smooth vector field v. This equation is said 
to be Hamiltonian if the vector field can be written as 

V (x) = O (x) {dH (x)) (2) 

where H (x) is the Hamiltonian function and fl (x) is the Poisson bi- vector 
(skew-symmetric, contravariant two-tensor) subjected to the Jacobi identity 

[O (x) , n (x)] = (3) 

defined by the Schouten bracket. In coordinates, if di = d/dx'^ the Poisson bi- 
vector is V, (x) = Q^'^ (x) dj A dk with summation over repeated indices. Then 
the Jacobi identity reads 

fi'^^iO'^'l = (4) 

where [jkl] denotes the antisymmetrization over three indices. It follows that in 
three dimensions the Jacobi identity is a single scalar equation. One can exploit 
the vector calculus and the differential forms in three dimensions to have a 
more transparent understanding of Hamilton's equations as well as the Jacobi 
identity. Using the isomorphism 

Ji — ^ijk^"^ ^5j)^~l?2,3 (5) 
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between skew-symmetric matrices and (pseudo)-vectors defined by the com- 
pletely antisymmetric Levi-Civita tensor e^fe we can write the Hamilton's equa- 
tions (2) in vector form 

v = JxVH (6) 
and in this notation the Jacobi identity (5) becomes 

J . (V X J) = (7) 

In this form, the Jacobi identity is recognized to be equivalent to the Probenius 

integrability condition for the vector J, or equivalently, the condition for the one 
form J = Jidx* to define a foliation of codimension one in three dimensional 
space [20], [21], [11]. 

A distinguished property of Poisson structures in three dimensions is the 
invariance of the Jacobi identity under the multiplication of the Poisson vector 
J (x) by an arbitrary but non-zero factor. More precisely, one can easily show 
that under the transformation 

J(x)-/(x)J(x) (8) 
of Poisson vector the Jacobi identity transforms as 

J • (V X J) {!{yi)f J • (V X J) (9) 
which manifests the invariance property. The identities 



J-v = 0, Vif-v = (10) 

follows directly from the Hamilton's equations (6), the second of which is the 
expression for the conservation of the Hamiltonian function. 

A three dimensional vector v (x) is said to be bi-Hamiltonian if there exist 
two different compatible Hamiltonian structures. In the notation of equation 
(6), this implies 

V = Ji X VH2 = J2 X Vffi (11) 

for the dynamical equations. The compatibility condition for Ji and J2 is 

defined by the Jacobi identity for the Poisson vector Ji + CJ2 for arbitrary 
constant c. Namely, Ji and J2 are compatible Poisson vectors provided they 
satisfy 

Ji ■ (V X J2) + J2 • (V X Ji) = 0. (12) 

The invariance properties of the Jacobi identity and the Hamiltonian functions 
enable one to extend the constant c to be a function of the conserved Hamiltoni- 
ans. More precisely, the Jacobi identity for the combination Ji + CJ2 of Poisson 
vectors gives 

(Ji X J2) • Vc = (Ji • (V X J2) + J2 • (V X Ji))c (13) 

which reduces to equation (12) whenever c is a constant. This linear equation 
can always be solvable for the function c resulting in considerable relaxation in 
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the compatibiUty condition. That means, locally every pair of Poisson vectors 
can be made compatible. 

It follows from the bi-Hamiltonian equations (11) and the identities (10) 
that Ji X V-ffi = J2 X VH2 = 0. That is, the Hamiltonian of one structure is 
the Casimir function of the other. Thus, a three dimensional dynamical system 
can be defined to be integrable if it is a Hamiltonian system with one Casimir. 
In this case, the flow can be represented by the intersection of surfaces defined 
by constant values of the integrals of motion (see [11] and references below for 
further details and examples). 

3 Jacobi Identity in Serret-Frenet Frame 

Let (t, n, b) denote the Serret-Frenet frame associated with a differentiable curve 
t x(t) in some domain of the three dimensional space R^. Throughout, 
V = {dx,dy,dz) will denote the usual gradient operator in local Cartesian co- 
ordinates. Given a vector field v, the unit tangent vector t, the unit normal n, 
and the unit bi-normal b can be constructed immediately as 

t(x) = ^ nW = p^{|^ b(x)=t(x)xn(x) (14) 

and they form a right-handed orthonormal frame except those vector fields v 
satisfying the condition imposed by 

t X (V X t) = 0. (15) 

It can be deduced from the vector identity 2t x (V x t) = V(t • t) + 2t • Vt = 

2t • Vt that this condition excludes essentially the flows with constant unit 
tangent and the points x at which the unit normal n (hence the bi-normal b) 
have zeros. That is, the cases one cannot have a Serret-Frenet frame. To avoid 
this we may assume that 

(Vxt)7^A(x)t (16) 

for arbitrary nonzero function A (x). That is, we exclude the dy- 
namical systems whose unit tangent vectors are the eigenvectors of the curl 
operator [22], [23]. 

We introduce the directional derivatives along the triad (t,n, b) as 

ds=t-V a„ = n • V ^6 = b • V (17) 

so that the variables (s,n, 6) are the coordinates associated with the Serret- 
Frenet frame. By inverting equations (17) we get the expression 

V = t^^ + ndn + hdb (18) 

for the Cartesian gradient in Serret-Frenet frame. Since t x V x t = t • Vt = dst 
the definition of the normal vector reduces to one of the Serret-Frenet equations 
justifying the name for the moving frame introduced [24]. 
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It follows from the identity in equation (10) that the Poisson vector J has 
no component along the unit tangent vector t. Hence, we set 

J = an + /3b (19) 

for unknown functions a (x) and /3 (x) satisfying +/3^ ^ 0. Using derivatives 
in Cartesian variables we find the expression 

J . (V X J) = (/3Va - aV/3) • t + a^n • (V x n) + (3% • (V x b) , . 

+ a/3 (n • (V X b) + b • (V X n)) ^ ' 

for the Jacobi identity. Assuming a and defining the function ^ = (j/a the 
Jacobi identity for a(n + /ib) gives 

t- V/u = n-V X n + /X (n-V x b + b- V x n) + /x^b-V x b (21) 

which is an equation involving only the unknown function /i. Obviously, this 
simplification is a manifestation of the invariance of the Jacobi identity under 
the multiplication of J by an arbitrary but non-zero function. Similarly, we 
may assume /3 7^ and define 77 = — l//x = — a//3 for which the Jacobi identity 
for the combination /3(b — r/n) becomes 

t-V?7 = b-V X b - (n-V x b + b-V x n) + jy^n-V x n. (22) 



To this end, we define the scalar quantities measuring the non-integrability 
(in the sense of Frobenius) of each unit vector in the Serret-Frenet triad 

fit = t- (V X t) , f7n = n- (V x n) , f^b = b- (V x b) (23) 

the first of which is necessarily not equal to any eigenvalue of the curl operator 
by the assumption 16. The integration of these quantities over the three space 
gives the so called helicities [25], [26] associated with the triad. We also introduce 
the sum 

n„b = n-V x b + b-V X n (24) 
of the cross-helicities for the normal and bi-normal vectors. 

Proposition 1 Let J = a(n + /ub) (or 3 = /3(b — ryn)^. Then, the Jacobi 
identity for J in moving coordinates is given by the Riccati equation 

dgH = On + M^nb + M^^b (25) 
(or dsT] = fib — T] flub + V^^a wUh ii=—l/r]) (26) 
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Thus, in the moving coordinates, the Jacobi identity or, equivalently, the 
existence of Poisson structure is expressible as a differential equation in arclength 
coordinates only. It may be interesting to note that the equation named after 
Jacopo Francesco Riccati originated from his investigations of curves whose 
radii of curvature depend only on a single variable [27] . The disappearence of 
the moving coordinates n and b from the Jacobi identity will become clear in 
the last section. In fact, they correspond to local conserved quantities and, as 
discussed in [11] for the globally integrable cases, may appear arbitrarily in the 
Poisson vectors. 

The Riccati equation (25) is equivalent to the linear second order equation 
2 / 9s fib 



d^u - i^-j^ + n„b ) dsU + On^bU = if fib 7^ (27) 
(or d^v - I - f^nb ) dsv + n^QbV = iiflny^O) (28) 



Q 

with the solutions being related by 

M=-^ ifJ^b^O (or r? = -%^ if 0„ ^ 0). (29) 

"b "n 

For the Poisson vectors of the above proposition at least one of the equations 
in (27) possesses two linearly independent solutions. 

The emergence of the Riccati equation as the Jacobi identity may be inter- 
preted as a relation between nonlinearity and superposition. The Jacobi identity 
as well as the compatibility condition for Poisson vectors are nonlinear restric- 
tions on some linear combinations of basis vectors. Both arc scalar equations in 
three dimensions. The Riccati equation, on the other hand, is known to be the 
only scalar equation admitting a nonlinear superposition principle [28]. 



4 Compatibility conditions 



To construct the bi-Hamiltonian structure, the Poisson vectors constructed from 
the linearly independent solutions of the Riccati equation must be compatible. 
That means, their linear combinations must also satisfy the Jacobi identity. 
Although, the multiplicative factors are left arbitrary in the construction of 
Poisson vectors they become important in the compatibility condition. Apart 
from the general result discussed in section two, wc shall restrict ourselves to 
the case where c is a constant in the combination Ji + CJ2. We shall show that 
the compatibility follows from the Hamilton's equations. 

First, we have the following result obtained by direct computation from the 
compatibility condition 

Proposition 2 Letai and ji^ he nonzero and be different functions for i = 1,2. 
For fib 7^ 0, the Poisson vectors J, = Q:j(n+/u^b) are compatible if 

ds\n — = (a*! - M2)^b- (30) 
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In obtaining Eq 30 we used the Riccati equations to eliminate derivatives of 
functions /U^ and /Z2. Next result shows that above equation is always satisfied, 
via Hamilton's equations, by Poisson vectors constructed from the solutions of 
the Riccati equation. 

Proposition 3 Let J = a(n + fih) and H define a Poisson structure for the 
dynamical system associated with v. Then 

a,lnU^-n-Vxb = /iOb (31) 
a 

Proof. With the Poisson vector in the assumptions we write the dynamical 
system as t =\\ v ||~^J x VH =|| v ||~^a(n+/ib) x V-ff. Taking cross products 
with n and b we get 

b = II V \\~\aVH - J(n • VH)) n = - || v \\~\anVH - J(b • VH)) 

(32) 

from which we obtain VH = a-^(l + /x^)--^(J(J-V-H') + J-"- || v \\).Here, we 
define J-*- = Q;(b— /in) for convenience. The integrability condition V x \/H = 
for the Hamiltonian function results, after taking dot product with J, in 



J ■ V X 



9slr^T^7r^ = TT7r—^ (33) 



where we used J • J-*" =0. The manipulations leading to the result is now 
straightforward and requires only the use of Jacobi identity for the derivative of 
II. ■ 

The proof of compatibility becomes obvious once we write equation 31 for 
each Poisson vector and subtract them. Similar results may be obtained for the 
Poisson vectors in the form J = /3 (b — Tjn) by assuming fi^ ^ and using the 
Riccati equation for the variable rj. 

We shall analyse the case fib = 0. This is the condition for the unit vector 
b to satisfy the Jacobi identity. On the other hand, Ob = reduces the Riccati 
equation to a linear first order equation resulting in one linearly independent 
solution for the construction of a Poisson vector of the form J = Q;(n+/xb). 

Proposition 4 The compatibility condition for Poisson vectors J = a(n+/ib) 
and b is dga + flnh = 0. This is satisfied by the Hamilton's equations. 
Proof. The equation can easily be obtained from the compatibility condition. 
Equation 31 with Ob = holds for the Hamiltonian structure for J. For the 
Hamiltonian structure with the Poisson vector h, we have t =|| v |P^b x VH 
for some Harriiltonian function H. Taking cross-product with h, solving for 
V-ff and taking the dot product of the equation resulting from the integrability 
condition V x ViJ = 0, we find 9slnl|v||=— b-Vxn which yields the result. 



Thus, the Poisson vectors obtained from solutions of Riccati equations are 
always compatible. 
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5 Bi-Hamiltonian structure 



Wc shall combine the results of the previous two sections on construction of 
Poisson vectors and their compatibility to present the main result. This will 
include the remaining case where the Poisson vectors are defined by the normal 
and bi-normal unit vectors of the Serret-Frenet triad. In connection with this 
particular structure we shall first relate the present work to the existing exam- 
ples of bi-Hamiltonian dynamical systems in the literature and then present the 
local form of the Nambu mechanics. 

Proposition 5 Every three dimensional dynamical system possesses two com- 
patible Poisson vectors. 

Proof. If both fin and Ob are non-zero, then any of the two Riccati equations 

which produce the same result by construction, give two Poisson structures com- 
ing from the linearly independent solutions of the corresponding second order 
equation. If flh = (or Sin = 0) then the first (the second) Riccati equation 
becomes linear with one linearly independent solution. Note that the sum, flnb of 
the cross-helicities of the normal vectors involves as an integrating factor in the 
integration for this Poisson structure. The other Poisson structure is defined by 
the bi-normal vector h (or the normal vector n) since Ob — (or On — 0) is the 
Jacobi identity for this. The compatibility conditions for these cases are shown 
to be satisfied via Hamilton's equations in the previous section. If, we have both 
On=0 and Ob = then n and b satisfies Jacobi identity and they become Pois- 
son vectors we sought. The compatibility condition is 0„b = 0. This implies 
that the function /x (orrj) must be a non-zero constant in the vector J for which 
the Riccnti equation is now trivial. The Hamilton's equations with the Poisson 
vectors n and b implies fJgln || v ||= n- V x b and Sgln || v ||= — b- V x n, respec- 
tively. Eliminating the vector v, we obtain the compatibility condition Clnb = 0. 
■ 

The last case of the theorem gives clues to relate the present work to the 

existing bi-Hamiltonian dynamical systems in the literature, in particular, to 
the Nambu mechanics. Recall that the three dimensional dynamical systems 
admitting bi-Hamiltonian structure with {3i,H2) and {32, Hi) are of the form 
of equation 11 with Ji and J2 satisfying Frobenius integrability conditions as 
Jacobi identities. In case that these vectors are globally integrable they are 
related to the conserved Hamiltonians by Ji = '.p^WHi z = 1,2 for some arbi- 
trary non-zero functions 9?^. In this case, the dynamical system has the form 
V = ij^VHi X VH2 first studied by Nambu in [2] . All explicitly constructed bi- 
Hamiltonian systems in three dimension has this form [6], [7], [11], [30]. The 
flow lines coincide with the intersections of the surfaces defined by contant val- 
ues of the Hamiltonians. Any constant appearing in this bi-Hamiltonian picture 
can be taken as arbitrary functions of i?i and H2. 

Our aim is to find the local version of this bi-Hamiltonian representation in 
three dimensions. This we shall do by first proving the more general result that 
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the local existence of a Hamiltonian function of the proposed form is equivalent 
to the existence of a Poisson vector (see also [29]). 

Proposition 6 There exists non-zero functions H and (p with 3 =ipVH when- 
ever J = n + fih satisfies the Jacobi identity. 

Proof. The condition V x ViJ = together with the form of J gives V x 
(n + fih) = a X (n + fih) where a = V In y^l + Ji? /f. Taking dot products with 
the unit vectors of the Serret-Frenet triad we obtain three equaiions one of which 
is algebraic and the other two expressing the s and n derivatives of the function 
jjL. All three contains terms involving the vector a. Eliminating this term for 
the s derivative results in the Riccati equation 25 for the Poisson vector J. The 
remaining two equations determine the function and the dependence of fx on 
the variable n. ■ 

Thus, finding two linearly independent solutions of the Ricatti equation com- 
pletely determines the local bi-Hamiltonian structure. 

In the local picture of the present work the normal coorclinatc^s n and b 
represents the local conserved quantities and they appear arbitrarily in the 
Poisson structures. The normal vectors n and b defining the bi-Hamiltonian 
structure in the last case thus corresponds, in the globally integrable ease, to 
the gradients of Hamiltonian functions defining Poisson vectors. Therefore, the 
manifestly bi-Hamiltonian equation t = n x b corresponding to the last case of 
the above proposition is the local version of the Nambu structure. 
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